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ABSTRACT 

A fibrous rare earth selective emitter is approximated as an infinitely long cylinder. The spectral emittance, is 
obtained by solving the radiative transfer equations with appropriate boundary conditions and uniform temperature. 
For optical depths, K R = a^R, where is the extinction coefficient and R is the cylinder radius, greater than 1 the 
spectral emittance is nearly at its maximum value. There is an optimum cylinder radius, R op| , for maximum emitter 
efficiency, r| E . Values for R ( are strongly dependent on the number of emission bands of the material. The 
optimum radius decreases slowly with increasing emitter temperature, while the maximum efficiency and useful 
radiated power increase rapidly with increasing temperature. 


INTRODUCTION 

Fibrous rare earth oxide selective emitters for thermophotovoltaic (TPV) applications have been of research interest 
for several years. Nelson (l) began working with fibrous emitters in the 1980’s. In addition fibrous emitters are 
being developed at Quantum Group (2) and the Auburn Space Power Institute (3). This paper develops the emittance 
theory for a fibrous emitter by approximating the emitter as an infinitely long cylinder. Since the fibrous emitters 
consist of bundles of 1 to 10 Jim diameter fibers this theory does not include the effects of the reflectance that 
occurs when radiation leaves a fiber and enters an adjoining fiber. The whole bundle of fibers is being approximated 
as a continuous cylinder. If the fibers are closely packed and the reflectance at the interface between a fiber and the 
medium in the voids between fibers is small, then the error resulting from the approximation should be small. 

The spectral emittance of the cylinder is obtained by solving the radiative transfer equation with appropriate bound- 
ary conditions. Knowing the spectral emittance allows the emitter efficiency to be calculated. As an example, 
emitter efficiency is calculated for an erbium-holmium aluminum garnet and thulium aluminum garnet (Tm^Al^O^) 
which are being studied at NASA Lewis. 


TEMPERATURE OF EMITTING CYLINDER 


As pointed out earlier (4), temperature drop across a planar or film type emitter causes a major reduction in the 
spectral emittance in the emission band of a selective emitter. However, in most cases for a cylindrical emitter there 
will be a negligible temperature drop. This can be seen by considering the steady state energy equation for an 
infinite cylinder with no internal heat generation and thermal conductivity, k th , and where we assume the tempera- 
ture, T, and radiation flux, Q, depend only on the radial coordinate, r. 


r 


dT 

k th — -Q 

dr 


= constant 


( 1 ) 


In order to avoid the term in brackets being singular at r = 0 it must vanish for all r. Thus at all r the conduction and 
radiation fluxes balance. 


^ih 


dT 

dr 


= Q 


( 2 ) 
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In other words, for all steady state conditions all the thermal energy bemg conducted into the cylinder at the outer 
radius. r= R. leaves the cylinder as radiation. 

The radiation flux, Q. will always be less than the blackbody flux o sb T', where a sb = 5.67xl0' 12 w/cm 2 K 4 is the 
Stefan Boltzmann constant. Therefore, define the following dimensionless variables. 


Q = 


Q 

tfstJs 4 




Where T s is the temperature at r = R. In this case equation (2) becomes the following. 


( 3 ) 


u 1 — 

-- = YQ (4) 

dr 

Where y is the ratio of the radiation flux to the thermal conduction flux. 

v - gshjjL _ a sh^s R ( c v 

• ~ T “ i tot 

k t h — k th 

in R 

For the ceramic materials used in most selective emitters, k th > 0.01 w/cmK. Also, for TPV applications T s < 2000K. 
Therefore, y < 4.5R(cm). So that ifR < 0.1 cm it is a reasonable approximation to neglect the right hand side of 
equation (4) and obtain the result T = constant (T = T s ). If y « l for a film or planar emitter then a linear tempera- 
ture variation results (4) rather than a uniform temperature. For the cylinder emittance calculation that follows, a 
uniform temperature is assumed. 



nf = index of refraction of cylinder 
n 0 = index of refraction of surroundings 
eosO ds = dr 


FIGURE 1. Schematic of emitting cylinder. 
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SPECTRAL EMITTANCE OF CYLINDER 


The spectral emittance is defined as follows. 


ei = ^ (R) 

X e bs (^,T s ) 


( 6 ) 


Where q^(R) is the radiation flux leaving the cylinder at r = R and e bs (>*,T s ) is the blaekbody emissive power where 
X is wavelength and T s is the cylinder temperature. 


e bs(^’T s ) = ni bs (A.,T s ) 


27thc0 

X 5 [exp(hc 0 /A,kT s ) - 1 


(7) 


Appearing in equation (7) is the blaekbody intensity i bs (^,T s ), w/cm 2 nm steradian, Planck’s constant, h, 
Boltzmann’s constant, k. and the vacuum speed of light, c Q . The radiation flux. q^(R), is obtained by solving the 
radiation transfer equations for the intensities, i^ + (R), and i^~(R). Where i^ + (R) is the intensity moving in the +R 
direction and i^~(R) is the intensity moving in the -R direction as shown in figure 1 . Assuming the intensities de- 
pend only on the radial coordinate leads to transport equations for i^ + (R) and i^~(R) identical to the planar case (5). 
These equations are written in terms of the optical depth, K, rather than the coordinate, r. 


K =«?l r 

(8a) 

* 

70 

II 

P 

>* 

70 

(8b) 


Where a^ is the extinction coefficient, assumed independent of r, and is the sum of the absorption coefficient, a^, 
and the scattering coefficient 0^. 


«A. =a ?i +0 X (9) 

The boundary conditions that must be applied are the following. At r = K = 0, from symmetry conditions. 

i+(0) = i^(0) at K = 0 (10a) 

At r = R or K = K R the intensity moving in the -R direction is equal to the reflected intensity. 

i^(K R ) = p to i+(K R ) at K = K R (10b) 

Where p (o is the reflectance at the cylinder outer radius, R. At r = R total reflectance occurs for certain angles of 
incidence, 0. At an interface between a material with a index of refraction, and a material with index of refrac- 
tion, nj, where > n^ radiation moving from i into j with angles of incidence 0 > 0 M , where 0 M is given by Snell’s 
Law will be totally reflected. Since n f > n Q , for the cylinder-air interface we have the following result for the 
reflectance, p fo . 


Pfo ” 1 


for 0 > 0jvi where 


Rm s cos2 0 M 


\ 

no 

n f V 


2 


For the case where 8 < 0 M ( u > (j. M ) we approximate p, Q by the reflectance for normal incidence (6). 


(11a) 
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Pfo = 


( 'Y- 

n f — n o 
n f + n o 


® < 0 M (P > Pm) 


(lib) 


Once the intensities i^ (K R .0) and i^ (K R ,0) are obtained the radiation flux, q^(K R ), can be calculated. 

«IX(K-r) = qJ(K R )-qx(K R ) 


( 12 ) 


t l^( K R) = 27t J Q=o i ^( K R- e ) cosesin9d0 = 27t J o i ^( K R.pH|i 

M^(Kr) = - 2 ^L i^(K R .e)cos0sin0d0 = 2n[ i^(0,n)ndn 
J v=nf l ^0 


(13) 


(14) 


Solution of the radiative transfer equations for q^(K R ) is presented in (5) for the film or planar case with a linear 
temperature variation through the film and for no scattering (a s = 0 in eq. (9)). Results for the uniform temperature 
cylinder can be obtained from these results by setting e fs = 0 (p fs = I) and AT = 0 in equations (33) to (36) of 
reference (5). When this is done the following result is obtained for the spectral emittance. 


n o0 Pfo) 

1 — 4E^(Kr) 


1-4E 3 (Kr) 

Pfo E 3( K R) + PM0 “Pfo) E 3 

Ml 

vPM J. 


uniform temperature and no scattering 


(15) 


Appearing in equation (15) is the exponential integral E^(x) defined as follows. 


E 


n 



n-2 


Z 



(16) 


The reflectance, p fo , is given by equation (lib) and |i M is given by equation (11a). 

As equations (15) and (11) indicate, for no scattering the spectral emittance of a uniform temperature cylinder 
depends on the optical depth, K R , and indices of refraction, n t > and n Q . For single crystal materials, such as rare earth 
doped yttrium aluminum garnet (YAG), scattering should be negligible. However, for polycrystalline rare earth 
oxides such as those being considered in references (!) to (3) scattering should be significant. In those cases 
equation (15) will overestimate the spectral emittance. 

Consider for the two limiting conditions K R = 0 and K R -4 The k R = 0 case corresponds to a wavelength 
where the material is transparent. While the K R —> °o case applies to an emission band of a selective emitter. Since 
E^(0) = 1/2 and x ^ 1 00 E(x) = 0, equation (15) yields the following: 


e*=0 for K R =0 

(17) 

lim ex = n 0 (l - p fo | 

(18) 

K R — >00 



Equation ( 1 8) is the usual result for an opaque body emitting (or absorl ing) in vacuum (n 0 = 1 ). 
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FIGURE 2 . Spectral emittance for cylinder of radius, R, at a constant temperature, T s , as 
a function of optical depth, K\ = ol\R, where is the extinction coefficient and nf is 
the cylinder index of refraction and n 0 is the surrounding index of refraction. 


In figure 2 the spectral emittance, e^, using equation (15) is shown as a function of optical depth, K R , for n 0 = 1 and 
n f = 1 , 1.5 and 2.0. As figure 2 shows increases rapidly with K R and reaches nearly its limiting value 
(eq. (18)) for K R = 1 . Notice also that for small K R as n f increases the spectral emittance rate of increase also 
increases. For most of the selective emitter materials, 1 .5 < < 2.0. 



2000 4000 6000 8000 

Wavelength, X nm 


lxUr 


FIGURE 3. Theoretical spectral emittance of Er<) 3 H02.7 AI 5 O12 cylinder of radius, 
R = 0.4 cm, calculated using measured extinction coefficient and index of refraction. 
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Once the extinction coefficient, a^, and index of refraction, n f , are known equation (15) can be used to calculate e^. 
Figure 3 shows the calculated spectral emittance for an erbium-holmium aluminum garnet cylinder (R = (3.4 cm) 
with 10% erbium and 90% holmium, Er Q ^Ho 2 7 Al s O p . This single crystal material is being considered for a film 
type selective emitter (5). The extinction coefficient and index of refraction were obtained from reference (5). 
Holmium has its main emission band centered at X ~ 2000 nm with smaller bands centered at X ~ 1100, 890 and 
750 nm. Erbium has its main emission band centered at X ~ 1500 nm with secondary bands centered at X ~ 1000, 
800 and 640 nm. All of these bands show up as regions of large emittance in figure 3. In the region 2000 < X 
< 4500 nm, Er 0 jHo 2 7 AI s O, 2 ls near ty transparent (a^ — > 0) and thus c ( is small. The highly oscillatory result in 
this region results from numerical error in a ^ (5). For the region X>X C = 5000 nm E x becomes large again. This 
large results from vibrational modes of the crystal lattice and is a characteristic of most rare earth selective 
emitter materials ( 1 ). We call the wavelength, X c , the long wavelength cutoff. 

EMITTER EFFICIENCY 

As a measure of the effectiveness of a selective emitter define the emitter efficiency as follows. 


*1E = 


X A, 

_ useful radiated power _ Qb _ Jo _ J 0 £ ^ e bs^T s )dX 

J 'OO <*oo 

q X (R)d\ e^e bs (X,T s )d?i 

o Jo 


total radiated power Qj 


(19) 


The numerator. Q b , is the power radiated in the wavelength region 0 < X < X In a TPV system X g corresponds to 
the bandgap energy, E g = hcJX^ of the PV cell. The denominator is the total radiated power, Q T . 



E 

o 

_c 

o 

Lm 

rj 

* 

O 

Cl 

*3 

o 


FIGURE 4. Emitter efficiency, iq E , and useful power, Q^, as ,t function of cylinder radius, 
R, at emitter temperature, T s = 1700 K for Erg 3 H02.7 AI5 ()|2- 


Consider how r| E will be behave as a function of the cylinder radius, R. As figure 2 shows, increases rapidly with 
optical depth, K R = a^R. Therefore, for the emission bands and the long wavelength region (large oc^ and 
X < Xg or X > X c ) will quickly approach its limiting value as R increases from zero. However, for the regions be- 
tween the emission bands and the long wavelength cutoff (small and X^<X< X c ) the spectral emittance will in- 
crease more slowly to its limiting value as R increases from zero. Thus, the numerator of equation (19), Q h , will 
rapidly rise from zero to its limiting value as R increases from zero. At the same time, the denominator, Q T , will 
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increase more slowly and will continue increasing with R while the numerator is increasing at a much lower rate. As 
a result, there will be an optimum radius, R opt , that will yield maximum Th* s ls i^ ustrale ^ > n figure 4 for 

Er 03 Ho-> 7 Al 5 O p where X g = 2200 nm (E g = 0.56eV) was chosen for use in equation (19) and T s = 1700 K. As can 
be seen, r| F rises rapidly to r) EMAX and then decreases slowly for R > R opr Also shown in figure 4 is the useful 
power radiated, Q b . As can be seen, Q b rises rapidly and then begins to level off. For T s = 1700K the optimum ra- 
dius is R opl = 0.34 cm. It should be mentioned that for radii of this magnitude the uniform temperature assumption 
becomes questionable. In figures 5(a) and (b) R op( , %(R opt ) and Q b (R opt ) are shown as functions of T s . Both 
n H (R op ,) and Q b (R op( ) * ncrease significantly with temperature while R opt decreases only 25 percent in going from 
T s = 1200 to 2000K. The large increase in x\ E results because the maximum value of the blackbody emissive power, 
e bs (^,T s ), shifts to shorter wavelengths as T s increases. Therefore, Q b increases faster than Q T as T s increases. The 
useful power, Q b , increases at least as T 4 ^ 



Emitter temperature, T s , K 



Emitter temperature, T s , K 

FIGURE 5. Optimum radius, R opt , maximum efficiency, tj EMax , and useful power, Q ft, at 
R = R opt as functions of emitter temperature, T s , for Erq 3 H02.7 AI5 Oi2- (a) Optimum 
radius, R 0 p t . (b) Maximum efficiency and useful power. 
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Now consider thulium aluminum garnet, Tm 3 AI 5 0 12 , as a cylindrical selective emitter. Thulium has its main emis- 
sion band centered at X~ 17()0 nm with smaller bands centered at X ~ 3 200, 800 and 700 nm. Thus compared to 
Er 0 3 H° 2 7 A1 5 0| 2 , which has 8 emission bands and thus a large region of high for X < Tm 3 Al s O p has a much 
smaller high emittance region for X < A, g . For Tm 3 Al s O j2 . X = 1900 nm was chosen. Just as for Er 0 3 Ho^ 7 Al 5 O p , 
thulium aluminum garnet has large emittance for X > X c = 5000 nm. 



Emitter temperature, T s , K 

FIGURE 6, Optimum radius, R opt , maximum efficiency, i) RMax , and useiul power, Qb> at 
R = R opt as functions of emitter temperature, T s , for Tm^ AI5 0\2- (a) Optimum 
radius, R 0 pt- (b) Maximum efficiency and useful power. 

In figure 6(a) and (b), R opt , r| R (R opl ) and Q b (R opt ) are shown as functions of T s for Tm 3 Al 5 O p . These results were 
calculated using measured values for a x and n f (5). The first thing to notice is that R op( for Tm 3 Al 5 O p is much 
smaller than R op{ for Er () 3 Ho 2 7 Al 3 O p . This occurs because thulium aluminum garnet has a smaller region of large 
emittance than erbium-holmium aluminum-garnet, as mentioned above Therefore, the majority of Q b results from 
emission in the main emission band centered at X ~ 1900 nm where oc^ = 30 cm -1 . Thus when R = 0.3 mm, K R = I 
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and will he near its maximum value (fig. 2). As a result, Q b will also be near its maximum and thus r| E will be a 
maximum. Therefore, we expect that R opt will occur near 0.3 mm for thulium aluminium garnet. Figure 6(a) 
substantiates this result. 

Comparing figures 5(b) and 6(b) we see that the emitter efficiencies are nearly the same for Er 0 ^H(v> 7 Al s O |2 and 
Tm^AI^Op. However, the useful power, Q b , is much larger for Er 07 Ho 7 7 A1 S 0 12 . This occurs for two reasons. First 
of all, erbium-holmium aluminum garnet has a much larger region of large emittance for X < X^ than thulium 
aluminum garnet. Secondly, X z = 1900 nm for thulium aluminum garnet whereas X^ = 2200 nm for 
erbium-holmium aluminum garnet. Thus more of the spectrum is included in Q b for Er^Hi^AhjO^- 

An important point about selective emitters is brought out by comparing figures 5(b) and 6(b). That is, r| E and useful 
power, Q b , do not behave in the same manner. Increasing Q b by adding more emission bands, as in the case of 
Er 0 7 Ho 0 7 AI s Oj 2 , does not mean that r| F will increase. 


CONCLUSION 

Most fibrous rare earth selective emitters consist of bundles of 1 to 10 fim diameter bundles. In this study that bundle 
has been approximated as an infinite cylinder. From the solution to the radiative transfer equations the 
spectral emittance of the cylindrical emitter was calculated. Several conclusions can be made about the cylindrical 
rare earth selective emitters. 

1 . For most rare earth selective emitters the temperature is uniform through the cylinder. 

2. When the optical depth K R = a^R > 1 the spectral emittance is nearly a maximum. 

3. There is an optimum value for the radius, R t , which yields maximum emitter efficiency. 

4. R opt strongly depends on the emitter material. For an emitter with only a single strong emission band R opt is the 

order of 0. 1 mm, whereas for an emitter with many emission bands R opt is the order of 1 mm. 

5. R op{ decreases slowly with increasing emitter temperature. 

6. The maximum efficiency. t1 EM aX’ anc * use ^ u * power, Q b , increase significantly with temperature. 
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